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UNDER  COMPRESSIVE  LOAD  (U) 

BY 

C  S  SMITH 


Summary  (U) 

A  review  is  made  of  the  problem  of  buckling  failure  in  fibre-reinforced 
plastic  panels  having  stiffeners  of  "hat"  section  under  compressive  load  acting 
across  the  line  of  the  stiffeners.  This  problem  is  of  practical  importance  in 
the  design  of  FRP  ships.  Alternative  analysis  methods  are  discussed,  including 
use  of  approximate  data  curves  and  folded-plate  analysis  to  determine  initial 
buckling  behaviour,  together  with  nonlinear  finite  element  analysis  referring 
to  imperfection  effects  and  post-buckling  behaviour.  Some  previously  unreported 
experimental  results  are  described.  These,  together  with  previously  published 
test  data,  are  correlated  with  theoretical  analysis  as  a  basis  for  making  up¬ 
dated  design  recommendations. 
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ABSTRACT  CARDS 


DESIGN  OF  TRANSVERSELY  STIFFENED  FRP  PANELS 
UNDER  COMPRESSIVE  LOAD  (U) 
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INTRODUCTION 


Fibre-reinforced  plastics  are  now  used  very  extensively  in  the  boatbuilding 
industry,  accounting  for  construction  of  about  80%  of  small  boats  up  to  25  m  in 
length  including  pleasure  craft,  workboats  and  fishing  boats.  FRP  construction, 
based  usually  on  the  use  of  polyester  resin  reinforced  by  E-glass  fibres  in  the 
form  of  chopped- strand  mat  or  woven-rovings  together  with  a  limited  amount  of 
unidirectional  tape  in  stiffeners,  is  also  widely  employed  in  larger  ships 
including  military  patrol  craft  up  to  40  m  in  length  and  minesweepers  of  up  to 
60  m  (1)*. 

2.  The  need  for  general  robustness  and  high  strength  under  transverse  load 
leads  in  many  cases  to  adoption  of  transversely  framed  hulls,  as  illustrated  in 
Figure  1,  containing  few,  if  any,  longitudinal  stiffeners.  Shell  and  deck  panels 
designed  in  this  way  are  clearly  less  efficient  than  longitudinally  stiffened 
panels  in  resisting  longitudinal  compression  associated  with  hull  bending: 
except  in  weight-critical  hulls  such  as  fast  patrol  craft,  however,  the  weight 
penalty  associated  with  transverse  framing  appears  to  be  preferable  to  the 
additional  complexity  and  hence  cost  of  fabricating  a  system  of  orthogonal 
stiffening.  Hull  designers  are  therefore  frequently  faced  with  the  need  to 
evaluate  the  compressive  strength  of  transversely  stiffened  deck  and  shell 
panels,  as  illustrated  in  Figure  2. 

3.  For  structural  design  purposes  a  ship's  hull  can  usually  be  treated  as  an 
assembly  of  stiffened  panels  which  may  be  examined  individually  assuming  either 
simple  support  or  appropriate  conditions  of  elastic  support  at  the  boundaries. 
Unlike  longitudinally  framed  panels,  whose  elastic  stability  is  inherently 
greater  and  in  which  a  substantial  post-buckling  reserve  of  strength  may  exist, 
transversely  stiffened  structures  are  likely  to  fail  catastrophically  at  loads 
approximately  equal  to  theoretical  critical  values  and  must  therefore  be 
designed  with  extreme  care  against  elastic  buckling.  Glass-reinforced  plastic 
(GRP)  panels  are  particularly  susceptible  to  instability  because  of  their  low 
material  stiffness  (Young's  modulus/ccxnpressive  strength  typically  <  80). 

Buckling  behaviour  is  complicated  by  the  anisotropic  nature  of  the  material  and 
the  "hat"  form  of  stiffeners  normally  used  in  GRP  hull  construction. 

4.  A  theoretical  study  carried  out  previously  (2)  has  identified  various  forms 
of  instability  occurring  in  transversely  stiffened  FRP  panels  and  has  established 
approximate  data  curves  and  accurate  computer-based  analysis  methods  for  evalu¬ 
ation  of  initial  buckling  stresses  and  mode  shapes.  Compression  tests  on  small- 
scale  perspex  and  GRP  models  (3,4,5)  and  a  limited  number  of  large-scale  GRP 
panels  (6)  have  provided  some  confirmation  of  theoretical  predictions.  The  pur¬ 
pose  of  the  present  report  is  to  assemble  all  available  experimental  data,  in¬ 
cluding  some  previously  unreported  results  of  tests  on  large-scale  GRP  panels, 

to  extend  theoretical  analysis  by  consideration  of  nonlinear  buckling  behaviour 
and  imperfection  sensitivity,  to  correlate  test  data  with  alternative  theoretical 
estimates  of  buckling  strength  and  hence  to  draw  up  recommendations  for  improved 
design  of  transversely  stiffened  GRP  panels  under  compressive  load. 


*(  )  =  References  on  page  16 
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ANALYSIS  METHODS 


5.  Assuming  simple  support  at  the  longitudinal  edges  (y  =  0,  B),  initial 
buckling  loads  Nxcr  and  corresponding  mode  shapes  for  a  longitudinally  com¬ 
pressed,  transversely  stiffened  GRP  panel  may  be  evaluated  using  the  methods 
outlined  below. 

Folded-Plate  Analysis 

6.  A  stiffened  panel  as  shown  in  Figure  2  is  treated  as  an  assembly  of 
rectangular  plate  strips  corresponding  to  frame  webs  and  tables  and  to  elements 
of  shell  laminate  contained  within  and  lying  between  hat-section  stiffeners. 

Each  plate  strip  is  related  to  x  and  y  axes  parallel  to  its  width  b  and  length  B 
respectively  and  is  assumed  to  have  principal  directions  of  elasticity  parallel 
to  its  sides  and  hence  to  satisfy  the  orthotropic  plane  stress-strain  relation¬ 
ship 
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Reinforcing  fibres  are  assumed  to  be  stacked  symmetrically  about  the  mid¬ 
thickness  of  the  laminate  so  that  no  direct  coupling  occurs  between  membrane 
and  bending  deformations:  the  moment -curvature  relationship  for  the  strip  is 
then 
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rigidities  per  unit  width  in  x  and  y-directions  and  D  =  G  h  /12  is 

xy  xy 

the  twisting  rigidity.  These  conditions  of  elasticity  are  effectively  satisfied 
by  most  marine-type  GRP  laminates  based  on  chopped  strand,  woven-roving  and 
unidirectional  reinforcement. 


7.  Equilibrium  and  compatibility  considerations  lead  to  the  following  pair  of 
equations  in  the  in-plane  displacements  u  and  v  of  a  plate  strip  subjected  to  a 
uniform  destabilizing  force  per  unit  width  N  in  the  x-direction. 


Similarly  the  equation  governing  bending  displacements  w  is  found  to  be 
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where  H  =  2  D  +  u  D 
xy  xx 


The  conditions  of  simple  support  at  y  =  0,  B  (ie  u  =  3v/3y  =  0  and 
2  2 

w  =  3  w/3y  =  0  at  y  =  0,  B)  allow  solutions  of  equations  3,  4  and  5  to  be 
assumed  of  the  form 
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in  which  a  =  ir/B  and  u^,  v^  and  w^  are  functions  of  x  independent  of  y. 

Substitution  of  (6)  in  equations  3,  4  and  5  leads  to  a  set  of  ordinary  differ¬ 
ential  equations  in  u^ ,  v^  and  w^  whose  solutions,  depending  on  the  boundary 

conditions  at  x  =  0,  b,  are  readily  found.  By  adopting  as  boundary  conditions 
unit  displacements  u^,  v^,  w^  and  rotations  dw^/dx  at  the  edges  x  =  0,  b  and  by 

obstitution  of  the  resulting  expressions  for  u,  v  and  w  in  equations  1  and  2, 
a  matrix  of  stiffness  coefficients  may  be  established  for  each  plate  element. 
Combination  of  element  stiffness  coefficients  under  conditions  of  continuity 
between  plate  strips  yields  a  set  of  4n  homogeneous  equations  which  express  the 
equilibrium  of  sinusoidally  distributed  moments  and  forces  at  all  n  plate  inter¬ 
sections  or  node  lines  in  the  idealized  structure.  These  equations  correspond 
to  a  standard  eigenvalue  problem  which  may  be  written  in  matrix  form 


where  K  is  a  square  stiffness  matrix  and  6^  is  a  column  matrix  defining  a 

buckling  mode.  In  general  an  infinite  number  of  eigenvalues  and  eigenvectors 
exist,  each  satisfying  equation  7  and  corresponding  to  a  particular  buckling 
mode.  Buckling  loads  may  be  found  by  varying  systematically  the  value  of 

and  evaluating  the  determinants  of  corresponding  stiffness  matrices  K:  critical 

values  N  are  those  for  which  the  determinant  of  K  is  zero.  This  process  can 
xcr 

be  carried  out  automatically  on  a  computer  with  buckling  loads  calculated  by 
iteration  to  any  required  accuracy.  When  a  critical  value  of  has  been  found 

the  corresponding  normalized  mode  may  be  obtained  by  setting  one  of  the  nodal 
displacements  or  rotations  equal  to  unity  and  reducing  the  homogeneous  equations 
(7)  to  a  set  of  one  fewer  nonhomogeneous  equations  which  are  solved  in  the  usual 
way  by  matrix  division.  Intermediate  modal  displacements  and  stresses  may 
finally  be  computed  at  any  point  on  each  plate  to  obtain  a  complete  definition 
of  the  buckling  mode. 

8.  The  analysis  method  outlined  was  developed  from  a  previous  generalized 
treatment  of  folded  plates  under  distributed  and  concentrated  lateral  loads  (7). 
Buckling  under  compressive  load  applied  in  the  direction  of  stiffeners  or 

biaxial  compression  may  be  examined  by  adding  appropriate  terms  to  equations  3, 

4  and  5  and  natural  frequencies  and  modes  of  vibration  may  also  be  evaluated  by 
including  inertia  terms  in  the  plate  equations  (8).  Analysis  of  this  type  has 
been  developed  independently  by  Wittrick  and  Williams  (9,10),  allowing  addition¬ 
ally  for  the  destabilizing  effect  of  shear  stresses  in  plate  elements  and  for 
somewhat  more  general  conditions  of  material  anisotropy,  and  has  been  applied 
more  recently  by  Viswanathan  et  al  (11). 

9.  Application  of  folded-plate  analysis  to  the  idealized  problem  shown  in 
Figure  2  has  established  that  the  lowest  initial  buckling  stress  usually  corres¬ 
ponds  to  a  local,  interframe  mode  having  one  of  the  forms  shown  in  Figure  3(a) 
(2).  Exceptionally,  in  structures  with  a  large  width  B,  these  modes  may  be  pre¬ 
ceded  by  Type-4  instability  having  either  (or  any  linear  combination)  of  the 
forms  shown  in  Figure  3(b).  Folded-plate  analysis  provides  an  economical  and 
theoretically  exact  method  of  examining  such  instability.  In  a  structure  of 
"infinite"  length  analysis  may  be  confined  to  a  single  frame-space,  or  half 
frame-space,  bounded  by  planes  of  synmetry  and/or  antisymmetry  depending  on  the 
mode  shape  under  consideration;  analysis  may  also  be  applied  to  a  structure  of 
finite  length  L  as  shown  in  Figure  2,  with  use  of  substructuring  to  achieve 
economy  of  computation  for  panels  containing  many  identical  frames. 

Data  Curves  for  Approximate  Evaluation  of  Interframe  Buckling 

10.  Guided  by  the  exact  buckling  displacements  and  associated  modal  stresses 
computed  using  folded-plate  analysis  approximate  models  of  interframe  buckling 
have  been  established  from  which  data  curves  may  be  derived  as  shown  in  Figure  4 
and  outlined  below. 

11.  Type-1  buckling  involves  bending  of  the  shell  laminate  restrained  by 
flexural  stiffness  of  the  frame  webs  and  tables  with  virtually  zero  displace¬ 
ment  at  intersections  between  the  shell  and  frames.  Provided  that  the  trans¬ 
verse  span  B  is  large  relative  to  the  frame  spacing  (b^  +  b^ ) ,  bending  in  the 

y-direction  may  be  ignored  and  the  shell  treated  as  a  continuous  beam  rotation- 
ally  restrained  at  points  of  attachment  to  the  frame.  The  buckling  load  per 


unit  width  N  is  then  independent  of  B,  depending  only  on  the  flexural 
rigidity  of  the  shell  in  the  x-direction,  the  dimensions  b  and  b^  and  the 
restraining  stiffness  k^  of  the  frame,  given  by 


k 


FI 


a(l 


where  a  =  4D  /b  and  c  =  a  +  2D  ,/b_.  Values  of  N  ,  obtained  by  solution 
x4  4  x3  3  xcr 

of  a  simple  eigenvalue  equation  referring  to  rotation  at  the  intersection  of  the 
shell  and  frame  web,  have  been  computed  for  a  range  of  k^  and  b^/b2  and  are 

plotted  non-dimens ionally  in  Figure  4(a)  from  which  Type-1  buckling  loads  may 
be  found  for  any  practical  frame-shell  configuration.  Assuming  b  >  b  it  is 

2  2  12 
evident  that  N  is  always  greater  than  w  D  , /b, . 

xcr  J  °  xl  1 


12.  Type-2  interframe  buckling  involves  bending  of  the  shell  restrained  by 
torsional  stiffness  of  the  frame,  with  a  modal  half -wavelength  equal  to  the 
frame  spacing  and  a  line  of  zero  displacement  at  the  centre  of  each  frame. 
Buckling  loads  again  depend  on  D^ ,  b^  and  b2  together  with  a  coefficient  k^2 , 

representing  torsional  stiffness  of  the  frame,  which  is  defined  by 


F2 


'  T7 

b2  B 


The  effective  torsional  rigidity  GJ  is  given  by 


Gj  =  4A*/*g  =  4A  \J 


b2  +  b3  2b4 
G  h_  G  H7  cTT 


where  Aq  is  the  area  enclosed  by  the  centrelines  of  the  frame  table  and  webs 

together  with  the  enclosed  portion  of  shell.  Evaluation  of  buckling  loads 

reduces  to  solution  of  an  eigenvalue  problem  with  two  degrees  of  freedom 

referring  to  displacements  and  rotations  at  frame/shell  intersections:  values 

of  N  ,  which  have  been  computed  for  a  range  of  k  .  and  b  /b_,  may  be  found 
xcr  F 2  1  2 

directly  from  Figure  4(b). 

13.  Type- 3  interframe  buckling,  which  involves  bending  of  the  shell  with  a 
medal  half-wavelength  equal  to  the  frame  spacing  and  lines  of  zero  displacement 
mid-way  between  frames,  is  restrained  by  bending  stiffness  of  the  frames  com¬ 
bined  with  a  membrane  shearing  action  in  the  shell  between  frames.  This  form  of 
instability  is  less  obviously  predictable  than  Types  1  and  2  and  is  not  a  form 
of  buckling  likely  to  be  encountered  in  panels  with  stiffeners  of  conventional, 
open  cross-section.  The  fact  that  the  preferred  mode  may  involve  flexural 
rather  than  torsional  deformation  of  the  frames  is  attributable  to  the  high 
torsional  rigidity  of  closed- sect ion  frames.  Buckling  loads  are  found  to  depend 
on  Dxl,  b^ ,  b2  and  a  stiffness  parameter  k^  which  represents  approximately  the 

restraining  influence  of  the  frame  and  is  given  by 


F3 


_  n_ 
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2  “V^l 

a  El  +  - - =■ 


4hlYl 


a  EAb1  -  41)^ 


where  is  the  shell  thickness,  EA  and  El  are  axial  and  flexural  rigidities  of 


a  composite  section  formed  by  the  frame  and  enclosed  portion  of  shell,  e  is 
the  distance  from  the  mid-plane  of  the  shell  to  the  centroid  of  that  section, 


a  =  tr/B  and 


G  + 


*2E  ,b? 

_ 

12 


Evaluation  of  buckling  loads  again  reduces  to 


solution  of  a  two-degree-of-freedom  eigenvalue  problem:  values  of  N  for  any 
k^  and  b^/b^  may  be  found  from  Figure  4(c).  xcr 


14.  Sufficient  information  is  given  above  for  use  of  the  data  curves  shown  in 

Figure  4:  a  more  complete  derivation  of  the  curves  may  be  found  in  Reference  2. 
Comparison  with  exact  folded-plate  solutions  for  a  range  of  geometries  has 
shown  (2)  that  these  curves  give  buckling  stresses  which  differ  usually  by  less 
than  10%  from  correct  values,  although  errors  may  exceptionally  be  as  much  as 
20%.  Type-2  and  Type-3  buckling  modes,  unlike  Type-1,  are  sensitive  to  the 
width  B  of  the  panel.  It  is  evident  from  Figure  4  and  should  be  noted  particu¬ 
larly  that  Type-2  and  Type-3  buckling  loads  N  may  be  less  than  the  value 

2  2  XCr 
ir  D^/b  :  the  assumption  that  interframe  buckling  corresponds  to  that  of  a 

simply  supported  strip  of  shell  of  width  b^,  which  has  previously  been  made  in 

design  of  transversely  stiffened  panels,  is  inadequate  and  unsafe. 

Finite  Element  Analysis 

15.  Finite  element  analysis,  employing  shell  elements  in  which  anisotropic 
material  properties  and  linear  destabilizing  (geometric  stiffness)  effects  are 
represented,  offers  an  alternative  means  of  evaluating  initial  buckling  stresses. 

Solutions  obtained  in  this  way  for  the  idealized  problem  of  Figure  2,  although 
less  accurate  and  more  expensive  in  computing  time,  have  been  found  to  corres¬ 
pond  quite  closely  with  folded-plate  solutions.  The  real  merit  of  finite 
element  analysis  lies  in  its  ability  to  deal  with  irregular  structural  geometry 
and  boundary  conditions:  buckling  solutions  extending  round  the  curved  portion 
of  structure  at  a  ship's  turn  of  bilge  as  shown  in  Figures  1  and  5  have  for 
example  been  obtained  in  order  to  check  boundary  conditions  assumed  in  the 
simpler  models  assumed  for  design  purposes.  As  described  below,  finite  element 
analysis  may  also  be  applied  incrementally  to  examine  the  nonlinear  buckling  and 
post-buckling  behaviour  and  imperfection  sensitivity  of  stiffened  panels. 

EXPERIMENTAL  RESULTS 


Previous  Test  Results 


16.  In  order  to  check  the  accuracy  of  analysis  methods  and  in  particular  to 
verify  theoretical  predictions  of  Type-3  instability,  compression  tests  were 
carried  out  on  two  small-scale  flat  rectangular  panels  (Models  PI  and  P2)  with 
transverse  hat-section  frames,  constructed  in  Perspex  (polymethylmethacrylate), 
as  illustrated  in  Figure  5(a).  In  order  to  check  design  assumptions  relating 
to  boundary  conditions  at  a  ship's  centreline  and  turn  of  bilge,  further  tests 
were  carried  out  on  small-scale  Perspex  models  of  a  vee-bottom  structure  incor¬ 
porating  a  stiff  keel  girder  (Model  P3  as  shown  in  Figure  5(b))  and  a  full  hull 
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compartment  (Model  P4,  Figure  5(c)),  Model  dimensions  (defined  as  shown  in 
Figure  6)  and  test  results  are  summarized  in  Tables  1  and  3.  In  Model  P3  the 
effective  panel  width  B  was  taken  as  the  distance  from  the  centreplane  to  each 
edge  while  in  Model  P4  B  was  taken  as  the  distance  fror  the  centreplane  to  the 
inside  edge  of  the  curved  bilge.  Mean  measured  properties  of  the  Perspex 

2  ? 

material  were  found  to  be  E  =  2850  N/mm  ,  y  =  0.37,  UTS  -  66  N/mm  .  A  full 
description  of  these  tests  is  given  in  Reference  3. 

17.  Small-scale  experiments  were  followed  by  compression  tests  on  three  large- 
scale  GRP  structures  (Panels  1,  2  and  3),  results  of  which  are  outlined  in 
Reference  6.  Tests  on  two  further  large-scale  GRP  models  (Panels  4  and  5), 
whose  results  have  not  previously  been  published,  are  described  below.  Tests 
have  also  been  carried  out  on  three  GRP  panels  (S2,  S3  and  S4)  representing  at 
about  1/4-scale  Panels  2,  3  and  4  (4,5).  Dimensions,  material  properties  and 
collapse  strengths  of  test  specimens  are  summarized  in  Tables  1,  2  and  3. 

Fanels  4  and  5 


18.  Test  Panels  4  and  5,  which  corresponded  to  the  bottom  shell  structure  i 
trial  ship  design,  were  nominally  identical  with  an  overall  length  of  609’  ~r 
and  effective  transverse  span  of  3048  mm.  Each  panel  was  stiffened  by  eight 
equally  spaced  transverse  frames  of  hat-section.  Fabrication  was  carried  . u' , 
following  standard  shipyard  procedures,  using  E-glass  reinforcement  and  a  r.  i.  ip  - 
type  isophthalic  polyester  resin.  The  shell  laminate  was  reinforced  by  3"  :  lie*- 

2 

of  800  g/m  woven  rovings  while  the  frames  incorporated  14  plies  of  800  •-  W! 

2 

with  12  plies  of  600  g/m  unidirectional  glass  in  the  tables. 

19.  Tests  were  carried  out  in  the  large  test  frame  at  ARE  in  a  rig  developed 
for  compression  testing  of  large  steel  grillages  (12).  As  shown  in  Figure  , 
panels  were  supported  horizontally  on  a  steel  platform  and  were  held  down  at 
their  ends  by  light  flexure  plates  and  along  their  sides  by  tie-bars  with  150  mm 
spacing.  This  form  of  attachment  was  designed  to  impose  conditions  of  simple 
support  at  the  ends  and  sides  of  test  panels,  with  ends  and  sides  free  to  move 
longitudinally  in-plane  and  sides  also  free  to  move  transversely  in-plane.  The 
ends  of  transverse  frames  were  restrained  against  torsional  rotation.  Longi¬ 
tudinal  compression  was  provided  by  seven  1500  kN  hydraulic  jacks,  evenly  spaced 
across  the  ends  of  test  panels,  acting  through  calibrated  hydraulic  load 
transducers.  Compression  was  applied  by  activating  the  jacks  at  one  end  of  the 
structure  from  a  common  pressure  source  and  reacting  loads  through  a  passive 
set  of  jacks  at  the  other  end.  In  order  to  distribute  jack  forces  into  test 
panels  in  a  uniform  manner,  avoiding  the  risk  of  premature  failure  close  to  the 
ends,  panels  were  reinforced  by  increasing  the  laminate  thickness  and  fitting 
steel  sandwich  plates,  bolted  and  bonded  to  the  GRP,  in  the  end  bays.  The 
sides  of  test  panels  were  also  reinforced  by  fitting  short,  discontinuous  steel 
sandwich  plates  between  transverse  frames.  The  ends  of  hat-section  frames  were 
strengthened  by  short,  built-in  steel  formers  through-bolted  to  external  steel 

p  iis. 


20.  Panel  5  was  tested  to  collapse  under  compressive  load  alone  while  Panel  4 
was  tested  under  a  combination  of  compressive  load  and  uniform  hydrostatic 
pressure.  Lateral  pressure  was  applied  by  means  of  a  water-filled  rubber  bag 
contained  between  test  panel  and  support  platform  as  shown  in  Figure  7. 

21.  Before  commencement  of  testing,  initial  deformations  of  each  panel  were 
surveyed:  maximum  interframe  displacements  w^  of  the  shell  laminate,  measured 

on  a  gaugelength  equal  to  the  frame  spacing,  together  with  maximum  overall 
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displacements  relative  to  the  ends  of  the  panels  measured  on  the  longitudinal 
centreplane,  are  listed  in  Table  1.  Local  distortions,  significant  in  relation 
to  interframe  buckling,  were  generally  very  small  (<  B/3800)  and  were  barely 
distinguishable  from  surface  irregularities  and  variations  in  laminate  thick¬ 
ness.  The  cov  (standard  deviation/mean  value)  of  laminate  thickness  t  was 

0.018  for  Panel  4  and  0.008  for  Panel  5.  During  tests,  displacements  were 
recorded  at  closely  spaced  positions  along  the  centreplane  of  each  panel. 

Strains  were  measured  on  both  sides  of  the  shell  laminate  and  on  frame  tables 
as  a  basis  for  vertical  alignment  of  jacks  and  as  a  check  on  the  uniformity  of 
compressive  load. 

22.  Panel  5  was  subjected  to  compressive  load,  applied  incrementally  with 
frequent  removal  of  load  to  allow  inspection  for  damage  and  to  provide  a  check 
on  permanent  set  and  loss  of  stiffness.  No  significant  damage,  permanent  dis¬ 
tortion  or  loss  of  stiffness  was  observed  up  to  an  average  compressive  stress 

2 

of  60  N/mm  .  At  this  stage  crackling  noises  commenced  and  as  load  was  increased 
whitening  of  the  shell  laminate  was  observed  at  the  base  of  frame  webs.  Large 
buckling  deformations  developed  as  shown  in  Figure  9  and  collapse  finally 

2 

occurred  explosively  at  an  average  stress  of  67.0  N/mm  taking  the  form  of  a 
fracture  over  the  full  width  of  the  panel  in  the  central  bay.  Failure  of  the 
shell  laminate  involved  delamination  with  shear  crimping  and  local  buckling  of 
the  delaminated  plies.  A  photograph  of  the  collapsed  panel  is  shown  in  Figure  8. 

o 

23.  Panel  4  was  sub}ected  first  to  a  uniform  lateral  pressure  of  103  kN/rT  , 
causing  overall  bending  deformation  with  an  amplitude  of  12.5  mm  together  with 
local  interframe  deformation  of  the  shell  laminate  of  maximum  amplitude  0.5  mm. 
Lateral  load  was  held  constant  at  this  level  while  compressive  load  was  applied 
incrementally,  simulating  a  design  condition  of  combined  compression  and  hydro¬ 
static  pressure.  No  damage  was  observed  at  compressive  stresses  of  up  to 

2 

60  N/mm  .  Collapse  finally  occurred  across  the  full  width  of  the  central  bay 

2 

at  an  average  stress  of  65.8  N/mm  .  Progressive  development  of  buckling 
deformations,  superimposed  on  displacement  caused  by  lateral  pressure,  is  shown 
in  Figure  10.  Failure  was  very  similar  to  that  of  Panel  5  and  lateral  pressure 
evidently  had  negligible  influence  on  compressive  strength. 

CORRELATION  WITH  THEORETICAL  BUCKLING  ANALYSIS 


24.  A  comparison  of  experimental  collapse  stresses  with  theoretical  initial 
buckling  stresses  for  Panels  4  and  5  and  for  the  other  GRP  and  perspex  panels 
is  made  in  Table  3.  Theoretical  buckling  stresses  were  estimated  in  three 
different  ways: 

(i)  from  the  data  curves  of  Figure  4; 

(ii)  by  application  of  folded-plate  analysis  to  a  single  interframe 
span  j 

(iii)  by  application  of  folded-plate  analysis  to  the  full  panel  length. 

It  is  evident  that  the  data  curves  provide  a  reasonable  first  approximation  to 
collapse  strength,  maximum  positive  and  negative  errors  being  +17%  and  -26%  of 
experimental  values.  Folded-plate  analysis  provides  a  more  accurate  estimate 
of  buckling  strength  with  maximum  errors  of  +5%  and  -20%.  In  the  latter  case 
(Panel  P3)  the  conservatism  of  the  theoretical  estimate  is  attributable  partly 
to  the  presence  of  a  torsionally  stiff  longitudinal  girder  at  one  edge  of  the 
panel  as  shown  in  Figure  5(b):  Panels  1,  3  and  S3  were  similarly  affected  by 
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the  presence  of  a  longitudinal  edge-girder.  Folded-plate  analysis  referring 
to  interframe  buckling  evidently  provides  a  close,  slightly  conservative 
approximation  to  the  lowest  initial  buckling  stress  for  a  complete  panel  of 
finite  length:  analysis  of  a  single  interframe  span  should  therefore  normally 
be  satisfactory  for  design  purposes. 

25.  It  may  be  noted  that  folded-plate  calculations  included  the  contribution 
of  frame  flanges,  which  effectively  increase  the  shell  thickness  over  a  span 
b  as  shown  in  Figure  6.  This  effect  was  omitted  (conservatively)  from  the 

data  curves  of  Figure  4.  Folded-plate  solutions  for  Panel  5  neglecting  frame 
flanges  gave  Type-1,  Type-2,  Type-3  and  Type-4  buckling  stresses  which  were 
respectively  2%,  1%,  8%  and  6%  less  than  values  including  flange  effects. 

26.  Correlation  of  theoretical  and  experimental  results  is  further  summarized 
in  the  form  of  scatter  diagrams,  shown  in  Figure  11(a)  and  11(b),  referring  to 
buckling  stresses  obtained  respectively  from  data  curves  (Figure  4)  and  from 
interframe  folded-plate  analysis. 

Modal  Analysis  of  Experimental  Buckling  Deformations 

21.  As  a  means  of  improving  the  correlation  of  theoretical  and  experimental 
buckling  behaviour  which,  as  illustrated  in  Figure  9,  may  be  complicated  by 
occurrence  of  several  closely  spaced  theoretical  buckling  stresses,  experi¬ 
mental  deformations  were  resolved  into  modal  components  using  the  orthogonal 

properties  of  theoretical  modes.  It  was  assumed  that  the  kt^1  experimental 
displacement  profile  w^  for  Panel  5,  measured  along  the  centreplane  as  shown 

in  Figure  9,  could  be  represented  by 


w  =  Z  a w.  .  (8) 

k  j=l  J 

where  a.  was  an  amplitude  coefficient  referring  to  the  j*'1  normalized 
theoretical  mode  .  The  orthogonality  condition  was  assumed  to  be 

L 

/  w .  w .  dx  =  0  (i^j)  .  (9) 

j  i  ,  ' 

o  J 

taking  x  as  the  distance  from  one  end  of  the  panel.  Coefficients  a  were 
then  found  from  the  expression 

L  lh  2 

a.  -  /  w .  w  dx  /  f  w.  dx  .  (10) 

]k  o  3  k  o  1 


integrations  being  performed  numerically  using  a  general  3-ordinate  Simpson's 
Rule.  Deflection  profiles  do  not  in  fact  fully  characterize  buckling  modes, 
which  involve  distortion  of  the  frames  and  in-plane  as  well  as  bending 
deformation  of  the  shell:  the  assumed  orthogonality  condition  is  therefore  an 
approximation.  It  was  however  found  from  a  numerical  check  on  orthogonality 
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distributed,  a  reasonable  permissible  level  of  N  for  design  purposes  might  be 

N  =  N  (g  -  3.1  a),  corresponding  to  a  probability  0.999  that  N  >  N  : 
xp  xcr  xu  xp 

dotted  lines  representing  these  load-levels  are  shown  in  Figure  11.  Because  of 

the  limited  number  and  scope  of  experimental  results,  however,  (influenced  by 

geometric  similarity  of  some  of  the  test  panels),  values  of  and  cannot 

be  regarded  as  sufficiently  accurate  to  provide  a  direct  basis  for  design  safety 
factors . 

36.  In  previous  ship  designs  (14,15)  a  factor  of  safety  of  4  or  more  has  been 
specified  against  calculated  elastic  buckling,  reflecting  uncertainty  about  the 
accuracy  of  calculations  and  the  relationship  of  initial  buckling  stresses  to 
structural  collapse.  Improved  understanding  of  buckling  failure  combined  with 
experience  gained  from  ships  in  service  indicates  that  this  requirement  is 
excessively  conservative.  As  a  compromise  between  past  practice  and  the 
approximate  statistical  criteria  discussed  above,  the  following  design  approach 
is  now  suggested :- 

(i)  where  the  lowest  initial  buckling  stress  Nxcr  is  estimated  accurately 

by  folded-plate  (or  equivalent  finite  element)  analysis,  a  per¬ 
missible  stress  N  =  N  /1.8  should  be  adopted:  this  factor  was 
xp  xcr  y 

chosen  on  an  "engineering  judgement"  basis  by  applying  a  factor  of 
1/1.5  to  a  load- level  at  which  significant  damage  might  be  expected 
to  commence,  judged  from  GRP  panel  tests  to  be  about  15%  below  Nxcr; 

(ii)  where  N  ^  is  estimated  by  use  of  the  data  curves  in  Figure  4,  a 

supplementary  factor  of  1/1.2  should  be  applied  to  account  for 
possible  inaccuracy  in  the  initial  buckling  stress,  giving  per¬ 
missible  N  =  N  /2.25. 
xp  xcr 

Recommended  permissible  stresses,  shown  as  full  lines  in  Figure  11,  evidently 
fall  below  mean  collapse  stresses  (corresponding  to  6.  and  $  )  by  5.1  a  and 
5.5  Oj  respectively.  12  1 

37.  It  is  assumed  in  the  foregoing  discussion  that  the  dominant  loading  on  a 
ship’s  hull  is  associated  with  wave  action  and  is  of  short  duration  (5  to  10 
seconds  at  most):  recommended  permissible  stresses  therefore  make  no  allowance 
for  creep.  Where,  exceptionally,  design  relates  to  a  dead  load  it  may  be 
necessary  to  introduce  an  additional  partial  safety  factor  accounting  for  creep 
effects. 

38.  It  should  be  noted  finally  that  theoretical  buckling  stresses  represented 
in  Table  3  and  Figure  11  and  reflected  in  the  above  discussion  of  design 
requirements  were  based  on  mean,  measured  elastic  moduli.  For  design  purposes 
a  supplementary  partial  safety  factor  accounting  for  variability  of  material 
properties  should  therefore  be  applied,  or  alternatively  lower-bound  material 
properties,  obtained  for  example  from  Reference  16,  should  be  adopted  in  design 
calculations. 
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32.  Results  are  summarized  in  Figure  14,  which  shows: 

(a)  compressive  load-shortening  curves  for  each  type  of  buckling; 

(b)  curves  of  average  compressive  stress  plotted  against  maximum 
lateral  displacement ; 

(c)  curves  of  maximum  outer-fibre  stress  in  the  shell  laminate  for  the 
lowest  (Type-3)  buckling  stress,  located  at  point  D  in  Figure  13(a); 

(d)  distribution  of  longitudinal  compressive  stress  across  the  shell 
laminate  at  various  load  levels. 

It  is  evident  that  each  type  of  buckling  involves  rapid  growth  of  lateral 
deformation  at  compressive  stresses  approaching  initial  critical  levels. 
Although,  as  in  elastic  columns,  a  stable  post-buckling  reserve  of  strength 
theoretically  exists,  such  behaviour  is  pre-empted  in  panels  of  practical 
geometry  by  material  failure,  resulting  from  high  bending  stresses,  before  a 
significant  increase  in  compressive  load  can  develop  beyond  the  initial  critical 
level.  Since  initial  distortions  of  moulded  FRP  structures  are  generally  very 
small  (those  of  Panels  4  and  5  were  somewhat  less  than  the  smaller  levels  of  w 

o 

represented  in  Figure  14),  large  lateral  displacements  and  hence  outer-fibre 
material  failure  do  not  occur  until  compressive  load  closely  approaches  the 
lowest  initial  buckling  stress:  this  critical  stress  therefore  provides  a  good 
estimate  of  collapse  strength. 

DESIGN  CONSIDERATIONS  AND  CONCLUSIONS 


33.  From  the  foregoing  discussion  it  is  evident  that  initial  buckling  stresses 
for  longitudinally  compressed,  transversely  stiffened  FRP  panels,  as  illustrated 
in  Figure  2,  may  be  estimated  accurately  by  folded-plate  analysis  referring  to 
local  inter frame  instability.  Care  should  be  taken  when  carrying  out  such 
analysis  to  examine  each  of  the  interframe  modes  shown  in  Figure  3.  A  reason¬ 
able  though  somewhat  less  accurate  estimate  of  initial  buckling  stress  may  be 
obtained  from  the  data  curves  shown  in  Figure  4.  Nonlinear  finite  element 
analysis  has  confirmed  experimental  observations  in  demonstrating  that  for  FRP 
panels  with  typical,  small  imperfections,  initial  buckling  stresses  obtained  by 
linear  analysis  provide  a  close  approximation  to  collapse  strength. 

34.  For  practical  design  purposes  it  is  suggested  that  preliminary  estimates 
of  initial  buckling  stress  should  be  made  using  the  data  curves:  at  a  later 
stage  in  design  these  should  be  checked  using  folded-plate  analysis.  Where 
folded-plate  analysis  is  not  available,  linear  finite  element  analysis,  although 
less  accurate  and  much  more  expensive,  may  be  used  instead  to  estimate  initial 
buckling  stresses.  Where  a  structure  is  highly  non-uniform  or  where  unusual 
boundary  conditions  exist,  finite  element  analysis  may  be  the  only  way  to  obtain 
accurate  results. 

35.  Specification  of  design  safety  factors  against  buckling  failure  should 
ideally  be  based  on  statistical  correlation  of  theoretical  buckling  stresses  and 
experimental  collapse  data.  Ultimate  compressive  strengths  N^u  may  be  related 

to  theoretical  initial  buckling  loads  Nxcr  by  an  empirically  determined  factor 

6  =  N  /N  .  From  the  theoretical  and  experimental  results  summarized  in 
xu  xcr  r 

Table  3  and  Figure  11,  mean  values  of  8  referring  to  data  curves  and  folded- 
plate  analysis  were  found  to  be  respectively  8^  =  1.079  and  §2  =  1.070,  with 

standard  deviations  =  0.125  and  a 2  =  0.093.  Assuming  that  8  is  normally 
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that  integrals  /  w.  w.  dx  (i  X  j)  were  usually  less  than  1%  and  in  almost  all 

°  1  3  L  j 

cases  less  than  10%  of  integrals  /  w.  dx. 

o  ^ 


28.  Assuming  in  accordance  with  linear  theory  that  modal  displacements  become 
infinitely  large  at  critical  stress  levels,  an  experimentally-based  estimate 
of  buckling  load  for  each  theoretical  mode  may  be  found  by  extrapolation  from 
a  plot  of  c  /a.,  against  the  experimental  compressive  stress  a  ,  as  illustrated 

X  JK  X 

in  Figure  12  for  Panel  5.  Theoretical  Modes  1  to  4,  shown  in  Figure  9,  were 
used  in  this  case  to  define  w..  From  Figure  12  it  is  clear  that  Mode  2 

J 

deformations  are  dominant  over  most  of  the  pre-collapse  load  range,  giving  an 
extrapolated  buckling  stress  which  is  consistent  with  the  theoretical  value  of 
2 

72.8  N/mm  .  Mode  4  deformations  are  also  present  and  the  extrapolated  buckling 

2 

stress  is  close  to  the  theoretical  value  of  90.3  N/mm  .  Mode  1  deformations 
are  not  evident  during  the  early  stages  of  load  application  but  emerge  strongly 
just  before  collapse.  No  meaningful  Mode  3  deformation  components  were 
identified. 


29.  The  procedure  outlined  above,  which  is  a  variant  of  the  Southwell  plot 
technique,  has  also  been  used  successfully  in  the  interpretation  of  buckling 
experiments  on  small-scale  perspex  models  (3). 

Nonlinear  Buckling  Analysis 

30.  In  order  to  obtain  an  improved  understanding  of  the  failure  process 
including  information  on  postbuckling  behaviour  and  imperfection  sensitivity, 
nonlinear  analysis  of  interframe  buckling  in  Panel  5  has  been  carried  out  using 
the  general  purpose  finite  element  program  ASAS-NL  (13),  implemented  on  the 
VAX  11/780  computer  at  ARE.  Three  different  finite  element  models  were 
examined,  as  shown  in  Figure  13,  in  order  to  represent  Type-3,  Type-2  and 
Type-1  buckling.  Eight-node  isoparametric  quadrilateral  shell  elements  were 
used  throughout.  Large  displacements  were  handled  using  an  updated  Langrangian 
(moving  co-ordinate)  formulation  with  application  of  a  modified  Newton-Raphson 
iterative  equilibrium  correction.  Anisotropic  material  properties,  as  defined 
in  Table  2,  were  assumed  to  remain  constant  at  all  load-levels. 

31.  Referring  to  Figure  13,  boundary  conditions  were  assumed  as  follows: 

(i)  plane  of  symmetry  at  AB  (u  =  3v/3x  =  3w/3x  =  0); 

(ii)  "moving"  plane  of  symmetry  at  CD  (3v/3x  =  3w/3x  =  0); 

(iii)  longitudinal  plane  of  symmetry  at  the  panel  centreplane  BD 
(v  =  3u/3y  =  3w/3y  =0); 

(iv)  simple  support  at  the  panel  edge  AC  (w  =  0,  otherwise 
unconstrained ) . 

Compressive  load  was  applied  incrementally  by  imposing  uniform  end-shortening 

displacements  6  on  the  plane  CD.  Initial  deformations  of  the  shell  laminate 

were  assumed  to  have  approximately  the  same  form  as  the  appropriate  interframe 

buckling  mode  with  sinusoidal  distribution  in  x  and  y-directions  and  an 

amplitude  w  :  in  each  case  two  distortion  amplitudes,  w  /B  =  0.0002,  0.001, 
o  o 

were  considered. 
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